Orthogonality catastrophe and shock waves in a non-equilibrium Fermi gas 
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A semiclassical wave-packet propagating in a dissipationless Fermi gas inevitably enters a "gra- 
dient catastrophe" regime, where an initially smooth front develops large gradients and undergoes 
a dramatic shock wave phenomenon. The non-linear effects in electronic transport are due to the 
curvature of the electronic spectrum at the Fermi surface. They can be probed by a sudden switch- 
ing of a local potential. In equilibrium, this process produces a large number of particle-hole pairs, 
a phenomenon closely related to the Orthogonality Catastrophe. We study a generalization of this 
phenomenon to the non-equilibrium regime and show how the Orthogonality Catastrophe cures the 
Gradient Catastrophe, by providing a dispersive regularization mechanism. 

PACS numbers: 73.22.Lp, 73.43.Jn, 78.70.Dm, 02.30.1k, 05.45.Yv 
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1. Introduction. When a Fermi gas is perturbed by 
a sudden switch of a local potential, it produces soft 
particle- hole pairs whose number grows as log(pj?i) with 
the size of the system. This phenomena is known as Or- 
thogonality Catastrophe 0] . It means that the overlap of 
the ground state of the Fermi gas with a localized poten- 
tial \B a ) (a state emerging as a result of a shake-off) with 
the unperturbed ground state (0| decays with the size of 
the system as (Q\B a ) ~ (ppL)~ a , where a = —S/tt and 
6 is a scattering phase of the potential. 

The effects of the Orthogonality Catastrophe are ob- 
served as the Fermi Edge Singularity - a power law reso- 
nance at the Fermi level occurring in transition rates in 
x-ray [2| or in tunneling experiments [2j. Recently this 
phenomenon has been exploited in a measuring device 
detecting local charge distribution in mesoscopic conduc- 
tors. 

The Orthogonality Catastrophe manifests itself differ- 
ently in systems out of equilibrium where energy relax- 
ation is small and electrons can diffuse out of the system 
without energy dissipation. The interest to tunneling out 
of equilibrium states is growing, but apart from recent 
works jj, |5( little is known. Perhaps the reason is that 
this problem can not be approached by methods tradi- 
tionally used for equilibrium states. 

2. Transition rates. Consider a non-equilibrium state 
(g | initially created in the Fermi gas, where we assume 
no interaction and ignore spin. This state evolves with 
the Hamiltonian 



2m 



(1) 



as (g{t)\ = (g\e~ ° . At time t we probe the state at the 
point a; by a sudden switch of a local potential U(x). 

In this letter we ask the following question. What is the 
probability to find the system in the ground state \B a {x)) 
of the the new, perturbed Hamiltonian H = Hq + U at 
time t. In other words, we are looking for a space-time 



dependence of the transition amplitude 
{g\e'^\B a {x)). 



(2) 



Such transition rates can be measured in transport ex- 
periments similar to those of Ref. 6]. There a quantum 
dot with a resonant level was brought into a proximity of 
a Fermi gas. When an electron tunnels out of the Fermi 
gas to the dot, the dot produces a potential seen by the 
Fermi gas a sudden shake-off. If the level suddenly be- 
comes unoccupied the transition amplitude reads 



(g\e-^^( x )\B' a (x)), 



(3) 



where \B' a {x)) is a ground state of the perturbed system 
with one extra particle. 

Other measurable quantities, e.g., tunneling current 
0,0] or generating functions of quantum noise , involving 
projections of evolving states (g(t)\ onto states of the 
perturbed gas other than the ground state, can also be 
computed using the methods developed below. 

Assume that the probing potential causes no backscat- 
tering and is well localized, so that a scattering phase S 
can be treated as a constant in the range of momenta of 
the wave-packet \g). Then a perturbed Hamiltonian is 

H = e "v(x) HQe -a V (x) 0| and itg ground state is 



\B a ) = e^\0) = (pfLY 



10}, 



S/tt, 



where the normal ordering separates the equilibrium 
part of the Orthogonality Catastrophe. Here (f(x) — 
2iri J x p{x')dx' is an antihermitian chiral Bose field and 
p(x) = ip^(x)ij}(x) is the fermionic density. Similarly, 
ip(x)\B' a (x)) - e^\B a (x)) - e^+^^lO), where we 
used a "bosonization" formula and set — 1 < a < 0. 

Summing up, we study space-time dependence of the 
transition rates 

r a = (g(t)\:e a ^:\0), r a+1 = (g(t)\ : e ^ a+1 ^ : |0), (4) 
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FIG. 1: [Color Online] A: A shock-wave solution of the Rie- 
mann equation 1101 . The dashed arrows indicate the velocity 
of the front. The vertical dashed lines are trailing and leading 
edges. The solid arrows show the relations between branches 
of the multi-valued solution of the Riemann equation and 
Whitham modulated particle P, hole Q and Fermi Pp mo- 
menta. B: Oscillations obtained by the Whitham method, 
the dashed red line shows the unphysical part of the Riemann 
solution. 

where we dropped space-time independent factors repre- 
senting the equilibrium part of Orthogonality Catastro- 
phe. We denote logarithmic derivatives as 

u = i—d x log u = —d x \ogr a T a+1 . (5) 

TO T a+1 TO 

We will see that the rates (@J) undergo complicated dy- 
namics, experiencing a shock-wave and a subsequent set 
of oscillations filling a growing spatial region. In fact, the 
shock wave occurs even at a = (or integer) , without the 
Orthogonality Catastrophe. However, its physics and the 
scale of oscillations are essentially different |8j . 

3. Semiclassical and coherent states. We will 
especially be interested in semiclassical wave- 
packets, i.e., states whose Wigner function 
W(x,p) = {g\e^ Px+Xp) \g), where P and X are 
momentum and coordinate operators, initially local- 
ized in the area of the phase space Ax 3> hpp 1 and 
Ap = \p — pf\ ^ Pf- This packet carries a large number 
of particles N — 2nhAxAp ^> 1. We may choose such a 
state to be coherent, i.e., given by {g\ = (0|e^p« A ^v^^ . 
This state corresponds to a smooth localized bump of 
electronic density as on Fig^and can be created by the 
action of a classical instrument. We also assume that the 
distance x between the initial origin of the wave-packet 
and a point of the measurement is large x 3> Ax. 

4- Hydrodynamic interpretation and the role of Orthog- 
onality Catastrophe. In the semiclassical approximation 
the amplitudes acquire a useful hydrodynamic inter- 
pretation. Let us assume that A PtP +k = Ak depends only 
on the momentum change, k, and write the initial state 
as (g\ = (olefv+frM*)^ where v+{x) = ^ k>Q A k Q t 
where V+(x) is an analytic function in the upper half- 
plane of x. In this case the density of the classical wave 
packet (g\p{x)\g) — — ^-lmVl_. On the other hand the 



initial values of the amplitudes are T a (x) = e aV +^ x \ and, 
therefore, initially ^{g\p(x)\g) = Reu(x). 

In the course of the evolution the above relation be- 
tween the density and the amplitudes is destroyed. How- 
ever, in the semiclassical limit and if a 7^ integer the rates 
P)l still contain all the hydrodynamic information. 

5. Dispersion of the electronic spectrum and non- 
linearity of the waves. It is commonly assumed that the 
linearization of the electronic spectrum at the Fermi sur- 
face H — Ef ~ 5Z p vf(p — PF^pipp captures the physics 
of the Orthogonality Catastrophe. If this were so, the 
time dependence of transition rates would be no different 
than its space dependence. The state (g(t)\ — {g\e~i VpPt 
simply translates the point of measurement: r(x,t) = 
t(x — Vf£) without any interesting dynamics. 

However, the approximation of linear spectrum is valid 
only for some time i « t c . It inevitably breaks down 
at larger time. The physics is simple: electrons in the 
denser part of the packet at the top of the bright side 
of the bump have higher momenta Sp = p — pf = hSp 
and, therefore, move with higher velocities v — vf = 
hhE" (pp)Sp than particles in front of them. Here 
^hE"(pF) =_J^ is a curvature of the spectrum at the 
Fermi point As a result, the wavefront steepens and 
eventually overturns (Fig. QJ. This is the shock wave we 
study using amplitudes Q as a "measurement". The 
results of the "measurement" depend on a and are espe- 
cially sensitive to whether a is an integer or not. 

The critical time of entering into the shock wave regime 
is about the time wave packet crosses the distance equal 
to the size of its front t c ~ We assume that t c is 

smaller than the ballistic time, so that dissipative effects 
in real systems do not have time to dissipate the shock. 

6. MKP equation of the soliton theory. Non-linear 
aspects of electron dynamics, can not be analyzed by el- 
ementary means. We have derived a fundamental equa- 
tion which determines both rates ijljl. It is the modified 
Kadomtsev-Petviashvili e qua tion (or MKP) - a known 
equation in soliton theory [lOj . Its bilinear form reads 

(<A - ^ D D T - ■ T -+i = °> ( 6 ) 

where D x f ■ g = f'g — fg' is the Hirota derivative. In 
fact, this equation holds for a more general class of matrix 
elements (g\e av ^\h), where \h) is any coherent state. We 
sketch the proof of the MKP at the end of this letter. 

Solutions of the MKP must be sought in the class of 
functions analytical in the upper half of the complex 
plane x. These are the properties of the matrix elements 
with respect to the Fermi vacuum - momenta of all ex- 
citations exceed the Fermi momentum. Analytical con- 
ditions are important. In particular they exclude soliton 
solutions of the MKP equation. 

In terms of |J5J we have another form of the MKP: 

it = ud x u + -^—d 2 x u. (7) 
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At a = 0, To = 1, and u — u. In this case a non- linear 
MKP equation becomes a linear Schrodinger equation 
for T\ = eT /" udx . Similarly, Q becomes a "complex 
Burgers equation" . Even in this simple case the dynamics 
of a semiclassical wave packet is not simple, but it can be 
obtained by elementary means studying the semiclassical 
limit of the solution of Shrodinger equation in the space 
of analytical functions in the upper half-plane. 

Analysis at a ^ requires methods of soliton theory. 

7. Multi-phase solution of the MKP equation. 10] . 
Assume that the initial state consists of a finite number of 
particles with momenta pi > pf and holes with momenta 
q l <p F , such that (g\ = (0\e^^ N . Then 



e"* detij(6ij + K a fa, qj)), where 



K a (pi,qj) 



sin(7ra) 



-4,, 



Pi - Pf 



a e i0i(x,t) 



pf - q% ) Pi - q 3 



(8) 



and0(p;,g;) = fa - qi)x - ^ (pf - qf)t, Bp = p F x-E F t. 
A formal solution for generic initial data is given by the 
determinant of the Fredholm operator 1 + K. 

This result can be obtained directly from the defi- 
nition of the matrix elements. One identifies the ker- 
nel K a (p,q) as the particle-hole amplitude K a (p,q) — 
(g\pq)(pq\e lPx - iHot \pq){pq\B a ), and A pq = (g\pq). The 
matrix element (pq\B a ) is the overlap between a particle- 
hole pair and the ground state of the perturbed Fermi gas 
computed in (pq\B a ) = (W)° for p jt q. 

Its singularity at the Fermi energy is a signature of the 
Orthogonality Catastrophe. 

In particular, the 1-phase solution is (p > pp > q) 



sin(7ra) / p — Pf\ °e" fi ' M ' 



PF 



(9) 



8. Quantum Riemann equation. In order to under- 
stand the MKP equation we recount the formulation of 
ID Fermi gas as quantum hydrodynamics, also known as 
(non-linear) bosonization, or collective field theory [Tlj . 
The quantum equation of motion of the chiral Fermi gas 
can be cast entirely in terms of the density operator 



d f u = ud x u, u = p. 

m 



(10) 



This is the quantum Riemann equation. This equation 
holds on coherent states generated by the density opera- 
tor. Its proof consists of a check that its l.h.s. commutes 
with all density modes pk = J e~ lkx p(x)dx. In their turn, 
the modes form the current algebra 

9. Classical Riemann equation and Shock waves. We 
note that the first two terms in the classical MKP equa- 
tion 10 are exactly the same as in the semiclassical ver- 
sion of the Riemann equation (|10J) . This is no surprise, 
since they are uniquely determined by the Galilean in- 
variance. One can neglect the third term in Q, or the 



quantum correction in (|10fl if the gradients are small. 
They, indeed, are assumed to be small initially. Also 
under a semiclassical condition one can neglect an imag- 
inary part of u. Then eq. p0|) becomes the Riemann 
equation of compressible hydrodynamics |l2| . 

Riemann's equation leads to shock waves: the velocity 
of a point with height u(x) is u(x) itself - higher parts 
of the front move with higher velocities. The bright side 
(ud x u < 0) of any smooth initial data gets steeper, and 
eventually achieves an infinite slope d x u(x Cl t c ) = oo at 
some finite time t = t c - a shock wave. 

After this moment the Riemann equation has at least 
three real solutions (Fig0 confined between X- (t) - the 
trailing edge, and x + (t) - the leading edge. They, and 
the critical point t c can be easily found from the im- 
plicit solution due to Riemann u(x, t) = f(x — u{x, t) ■ t), 
where f(x) = u(x,0) is the initial profile. For a typ- 
ical wave packet with a height Ap/m and width Ace 
a critical time is of the order of t c ~ mAx/Ap. The 
leading edge (in the Galilean frame moving with veloc- 
ity Dp) moves with velocity Ap/m: x+ ~ (Ap/m)t. If 
f(x) ~ x~ n at i > Ax the trailing edge delays, pro- 
gressing as x- ~ Axlt/tc) 1 /^ 1 " 1 . 

Let us label the solutions of the Riemann equation at 



,(0) 



is a single-valued solution outside the shock 



wave interval, u^>u^ >u^ are three ordered solutions 
in the shock wave interval x~ (t) < x < x + (t) . The branch 
smoothly merges with at the trailing edge, while 



smoothly merges with at the leading edge. 



10. Dispersive regularization and the role of Orthogo- 
nality Catastrophe. Obviously, the approximation lead- 
ing to the Riemann equation fails when gradients are 
large. Then, the neglected gradient terms become impor- 
tant. They regularize the "gradient catastrophe". The 
regularization and the subsequent physics of the shock 
wave is very different for an integer a = 0, and in the 
case of the Orthogonality Catastrophe, where a is irra- 
tional (a rational a involves some additional structures). 

If a = the solution of the Schrodinger equation at 
t > t c is well approximated by u' ' at x < X-(t) with 
an abrupt fall to at x > X-(t). In the case of the 
Orthogonality Catastrophe (a ^ integer), we face com- 
plexity of the nonlinear equation J7J). In this case the 
entire interval X-<x <x+ is filled by oscillations (Fig[IJ. 

11. Whitham modulation. Despite the integrability of 
the MKP equation, its solution in the form of the Fred- 
holm determinant, with the kernel JSJ, is rather compli- 
cated and the initial value problem is generically difficult 
to solve. Fortunately, a powerful approximate method to 
describe the shock waves has been developed in a semi- 
nal paper 0] . The method suggests to glue a solution of 
the Riemann equation which is valid for x < X- (t) and 
x > x+(t), to a periodic solution. In the first approxima- 
tion the 1-phasc solution can be used. The amplitude 
and the period of the wave have to be modulated in order 
to match very different values of the front at the trailing 
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x_ (t) and the leading x+ (t) edges of the shock. 

Modulated non-linear waves are the subject of the 
Whitham theory 01 . The latter states that modulated 
waves have the form of a multiphase solution ©, which 
moduli and phases are smooth functions of space-time 
p, q,p F , 6, 9 F -» P(x, t), Q(x, t),P F (x, t), Q(x, t), Q F (x, t). 

In our case the moduli have a clear physical interpre- 
tation. They are momenta of soft particle-hole pairs pro- 
duced at the Fermi level, which is also changing in space- 
time. The phases in (jHJ obey the Whitham equations: 

e=E(py-E(Q),a x e=p-Q, e F =E(p F ),d x e F =p F , 

where E(P) — £- is a modulated energy. The Whitham 
equations for the moduli are determine d by Galilean in- 
variance. They are again Riemann eqs. [15J 

P + d x E{P) = Q + d x E(Q) =P F + d x (P F ) = 0. (11) 

The initial data for the Whitham equations are chosen 
so that the 1-phase oscillatory solution © is glued to a 
(non-oscillatory) solution of the Riemann equation 
at the leading and the trailing edges. 

Assuming that A pq is smooth, we notice that the 1- 
phase solution stops to oscillates when a hole is absorbed 
at the Fermi level at the trailing edge, and when a par- 
ticle is created at the Fermi level at the leading edge 
Q(x-) = P F {x-), P(x + ) = P F (x + ). This yields that 
at the trailing edge u(x-) — (i/m)d x (Q + <d F ). Accord- 
ing to the Whitham equation id x (Q + Q F ) = P — Q + P F . 
which is just P at x = x- . Therefore uW = P are the 
boundary data for the Whitham equation for P at the 
trailing edge. Since vS 1 ^ is a solution of the Riemann 
equation P — holds in the entire oscillatory interval. 

At the leading edge, u — (i/m)d x (<d F — 6) — 
(l/m)(P F +Q-P) = Q/m. Therefore, u^(x + ) = Q(x + ) 
is a boundary data and also a solution of the Whitham 
equation for Q. In a similar fashion one concludes that 
the modulated Fermi momentum P F is given by the 
branch of the Riemann equation. Summing up 

x-(t)<x<x+(t) : P = u {1) < P F = u {2) < Q = u (3) . 

Being substituted into the 1-phase solution 10 these for- 
mulas give an explicit (approximate) solution of the MKP 
in the oscillatory region shown in Fig. ^ fl^t 

12. Derivation of the MKP equation. We sketch 
the derivation of the MKP for the amplitudes r a = 
(g\e a ^ x '^\h) where \h) is a generic coherent state (we 
focused on \h) = |0) in the paper). For a more detailed 
discussion of the relation between the dynamics of the 
Fermi gas and soliton theory see [l(ij . 

First, with the help the quantum Riemann equation 
(|10fl . we compute the action of the Schrodinger operator 

(idt + ^~d 2 x ) : e av := a(a + 1) : e aip T :, 
2m 

(-idt + -^-d 2 x ):e^ a+1 ^:= a(a + 1) :e^ +1 ^f:, (12) 



where T —: ip 12 : —if" and T —: ip' 2 : +cp" are holomor- 
phic (antiholomorphic) components of the stress-energy 
tensor of a chiral Bose field. 

Using these formulas we write the eq. (0 in the form 

(g\e a v\h) + {g\e^ a+1 ^\h) ~ (g\e a ^\h) {g\e( a+1 ^\h) 

where J = d x tp — 2-nip is the current of the Bose field, 
and the expressions are understood to be normal ordered. 

In terms of fermions T{x) ~: ip^ {x)d 2 ip(x) :, and using 
the bosonization formula : e av :~: e^ 1 ^^ {x)il){ioo) : 
we rewrite the numerator of the first term in the l.h.s. as 
a four fermion insertion: 

lim d 2 (g\ :^{z)i;Hx)^(y)^(ioo)e aip{x) : \h). 

z.y — >x y 

One may now apply Wick's theorem, carefully taking care 
of normal ordering, to write this in terms of matrix el- 
ements with two fermion insertions. Then, writing the 
fermions in terms of the Bose field and taking the limit 
one proves the MKP equation (0. 

13. Summary. The density of the semiclassical wave 
packet measured by a sudden switching of the poten- 
tial initially behaves according to the classical Riemann 
equation, i.e., similar to a propagating disturbance in a 
classical compressible liquid. At the time t c the wave 
packet enters a shock wave regime. It collapses emanat- 
ing particle- hole pairs resulting in modulated oscillations. 
The wave vector of oscillations is of the order of Ap - the 
"height" of the initial wave packet. It is much smaller 
than the Fermi scale p F . The oscillations occupy an in- 
terval whose leading edge propagates with the velocity 
exceeding the Fermi velocity by Ap/m. The oscillations 
are a distinct signature of the Orthogonality Catastro- 
phe. An observation of quantum shock waves, say, on 
the edge of Integer Quantum Hall Effect would be yet 
another manifestation of the quantum coherence. 
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